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i 1 .  ‘INTRODUCTION  AND  SUMMARY 
Consider  the  classical  MANOVA  model 


Also,  and  are  of  order  P^Pj  and  nxr..  respectively,  P  *  Pj  +  P2  and 
r  =  rl  +  r2‘  *n  Edition,  Y.  and  E.  are  of  order  nxp^. .  Then,  we  are  interested 
in  testing  the  following  hypotheses: 


Problem  [I]  H:e12  *  0  and  r12  a  0  vs  K:  not  H 

Problem  [II]  H:e12  =  0,  e2^  =  0  and  Z^  =  0  vs  K:  not  H 
Problem  [III]  H:Zi2  *  0  under  e12  =  0  and  021  =  0 

vs  K:Ij2  ?  0  under  e12  =  0  and  e21  =  0. 

The  motivation  behind  each  problem  is  stated  in  Section  2  and  some  examples 
are  also  given  there.  In  this  section,  some  formal  features  of  the  problems 
are  made  clear  and  our  results  are  briefly  summarized  together  with  some  results 
in  the  literature.  A  basic  feature  in  the  problems  treated  here  is  that  in  each 
hypothesis  the  independence  (z^2  =  0)  between  Y1  and  Y2  is  included  corresponding 
to  the  structure  of  the  regression  coefficient  matrix  0. 


K-: 
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In  Problem  [I],  the  hypothesis  will  be  regarded  in  Section  2  as  a 
formulation  of  the  hypothesis  of  no  causality  from  Xj  to  Y2  where  X  =  [Xj,X2] 
may  be  random  but  is  fixed  with  full  rank.  Also  e12  =  0  in  the  hypothesis 
may  be  viewed  as  a  special  case  of  the  GMANOVA  (general  MANOVA)  hypothesis 

M10M2  =  0  (1.3) 

where  M^  and  M2  are  fixed  matrices  of  full  rank.  In  fact,  Mj  =  [1,0]  and 
M2  =  (i)  imP^ies  ei2  =  The  Pro^em  of  testing  (1.3)  in  the  GMANOVA  model 
of  the  form  Y  =  Z^0Z2  +  E  is  known  as  the  GMANOVA  problem  and  has  been  treated 
by  Potthoff  and  Roy  (1964),  Rao  (1965,1966),  Khatri  (1967),  Krishnaiah  (1969), 

Gleser  and  01  kin  (1970),  Fujikoshi  (1973),  Kariya  (1978),  Marden  (1982a)  and 
others. Potthoff  and  Roy  (1964)  proposed  the  GMANOVA  model  and  considered  adhoc 
procedures  for  testing  the  general  linear  hypothesis  on  the  location  parameters. 

Rao  (1965,1966)  reduced  the  problem  to  testing  the  general  linear  hypothesis 
under  a  conditional  model.  Khatri  (1966)  derived  the  LRT  for  testing  the 
general  linear  hypothesis  under  Potthoff  and  Roy  model .  Later;*  Gleser  and  OTkiri.  ( 197( 
gave  a  canonical  reduction  of  the  problem  and  discussed  the  LRT  procedure  using 
the  canonical  form.  However  they  have  not  treated  Problem  [I].  It  should  be 
noted  that  when  (1.3)  is  dealt  with,  the  presence  of  a  general  M2  in  (1.3) 
affects  the  covariance  structure  so  that  in  the  case  of  (1.3)  should  be 

replaced  by  the  hypothesis 

Z  =  M2TM2  +  P2AP2  (1.4) 

I 

where  P2  is  a  matrix  of  full  rank  satisfying  P2M2  =  0  (see  Kariya  (1985),  pp 
175-176).  That  is.  Problem  [I]  is  considered  equivalent  to  the  problem  of  testing 
(1.3)  and  (1.4)  simultaneously  since  both  problems  give  the  same  canonical  form. 
Hence,  solving  the  former  problem  implies  solving  the  latter 
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problem  and  vice  versa.  Now  for  Problem  [I],  applying  the  invariance  principle, 
we  first  expand  the  power  function  of  an  invariant  test  in  the  neighborhood  of 
the  hypothesis  and  based  on  it  propose  a  test,  in  Section  3,  which  maximizes 
the  power  in  a  slightly  restricted  neighborhood  of  the  null  hypothesis.  Of 
course  due  to  the  local  optimality,  it  is  admissible.  The  test  statistic  there 
is  a  linear  combination  of  the  LBI  test  statistic  for  testing  =  0 
(Schwartz  (1967))  and  the  LBI  test  statistic  T£  for  testing  Z^  =  ®  under 
*12  *  The  latter  test  is  equivalent  to  the  LBI  test  of  independence  with 
some  data  missing  in  Eaton  and  Kariya  (1983).  Because  of  the  form  of  T^,  Tj 
and  Tg  are  correlated  and  so  our  test  is  not  equal  to  a  test  combining  the  two 
independent  LBI  test  statistics  Tj  and  T3  for  the  two  separate  hypotheses 
e^2  =  0  and  Z^  =  0  (without  6^  =  0),  though  T^  is  the  same.  This  is  a  feature 
of  the  joint  treatment  of  the  two  hypotheses  and  it  implies  that  a  test  combining 
the  two  independent  statistics  T,  and  T3  which  are  LBI  for  each  hypothesis  does 
not  maximize  the  local  power  in  any  direction  except  for  the  case  that  the  test 
depends  on  Tj  only  and  that  the  alternative  space  is  restricted  to  the  space  on 
which  Zjg  *  0-  The  problem  of  how  to  combine  independent  tests  is  discussed  in 
the  literature  (see  e.g.,  Marden  (1982b)),  though  we  do  not  discuss  it  here. 

But  the  LRT  statistic  for  Problem  [I]  gives  a  natural  combination  for  two  sepa¬ 
rate  hypotheses.  In  fact,  it  is  the  product  of  the  two  independent  LRT  statistics. 
This  might  support  the  idea  that  we  separately  treat  the  hypotheses  and  then 
combine  the  two  tests.  However,  as  has  been  observed  in  Eaton  and  Kariya  (1983), 
even  when  =  0,  the  LRT  for  testing  z^  =  0  ignores  the  additional  information 
(data)  available  through  =  0*  In  this  sense,  the  above  fact  may  not  be 
seriously  taken  into  account.  The  asymptotic  null  distributions  of  the  test 
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based  on  and  anc*  the  LRT  are  derived  in  Section  5  and  the  unbiasedness 
of  the  LRT  is  shown.  It  is  noted  that  the  group  leaving  the  problem  invariant 
is  small  so  that  the  power  function  of  an  invariant  test  including  the  LRT 
depends  on  many  parameters  including  the  canonical  correlations. 

In  Problem  [II],  the  hypothesis  will  be  regarded  in  Section  2  as  the 
hypothesis  of  no  additional  information  in  canonical  correlation  analysis  or 
a  formulation  of  the  hypothesis  of  no  causality  from  to  and  from  X2  to 
Yj,  where  X  may  be  random  but  is  fixed  with  full  rank.  Here  the  restrictions 
0j2  =  0  and  ~  0  are  special  cases  of 

MjGMg  =  0  and  M3eM4  =0.  (1.5) 

However,  the  two  GMANOVA  type  restrictions  in  (1.5)  cannot  be  expressed 
as  a  single  GMANOVA  hypothesis  of  the  form  M^eM2  =  0.  That  is,  the 
problem  of  testing.,  (1.5)  even  in  our  MANOVA  model  Y  =  Xe  +  E  is  no 
longer  the  GMANOVA  problem,  and  difficult  to  treat  unless  M^ 

i  iii 

and  M3  are  nested  relative  to  X  X  or  orthogonal  relative  to  X  X,  i.e.,MjX  Xm3  = 
0  (see  Kariya  (1985)  p  143).  Since  Mj  =  [1,0]  and  M3  =  [0,1]  in  our  present 
case,  MjX  XM3  =  X^X2  and  hence  without  XjX2  =  0  the  problem  of  testing. 

Problem  [IV]:  H:0j2  3  0  and  02j  =  0  is  difficult  to  treat.  In  fact,  it  is 
not  only  difficult  to  derive  the  LRT  explicity  but  it  is  also  difficult  to 
find  a  similar  test  detecting  both  0^  =  0  and  0^  =  0  in  a  meaningful  manner 
(see  Section  7).  On  the  other  hand,  the  hypothesis  on  the  covariance  structure 
which  corresponds  consistently  to  the  hypothesis  (1.5)  is  expressed  as 


1 2 r i^ 2  *  ^2^1^2  a^  ^  =  ^4*2^4  ^4^2^4' 


(1.6) 


where  P2M2  =  0  and  P^  =  0.  Since  in  Problem  [II]  M2M4  =  0,  we  can  take 
?2  =  and  P^  =  so  that  the  two  covariances  in  (1.6)  become  the  same. 

I 

In  Section  4,  first  in  the  case  of  X^X2  +  0  we  analyze  Problem  [II]  via 
invariance,  but  because  the  group  leaving  the  problem  invariant  is  quite 
small,  no  sufficient  reduction  is  obtained  and  the  space  of  a  maximal  invariant 

I 

parameter  is  of  high  dimmension.  Hence  in  the  case  of  X^X2  f  0  we  simply  show 
the  unbiasedness  of  the  LRT  derived  by  Fujikoshi  (1982).  The  LRT  statistic 
here  is  the  product  of  the  three  LRT  statistics  for  the  three  separate  hypotheses 
e  12  =  0,  =  0  an<*  si2  *  ®  but  the  three  are  dependent.  It  is  noted  that  it 

is  difficult  to  consider  the  monotonicity  of  the  power  function  of  the  LRT 

I 

because  of  the  high  dimmensional  parammeter  space.  Next  in  the  case  of  X^X2  =  0, 
we  expand  the  power  function  of  an  invariant  test  in  the  neighborhood  of  the 
null  hypothesis  and  based  on  it  propose  a  test,  which  is  a  linear  combination 
of  three  statistics  Rj,  R2  and  R3.  Here  similar  to  Problem  [I],  Rj  and  R2  are 
respectively  the  LBI  test  statistics  for  the  hypotheses  e12  =  0  and  e21  =  0 
and  R3  is  the  LBI  test  statistic  for  Z12  «  0  under  e12  *  0  and  021  =  0.  Hence 
Rj,  R2  and  R3  are  dependent.  This  is  a  feature  different  from  the  separate 
treatment  of  the  three  hypotheses.  The  asymptotic  null  distributions  of  this 
test  as  well  as  the  LRT  are  given  in  Section  5. 

Problem  [III]  was  treated  by  Kariya,  Fujikoshi  and  Krishnaiah  (1984) 
(abbreviated  as  KFK  henceforth).  In  this  model  X  is  fixed  but  it 
may  not  be  of  full  rank.  The  model  (1.1)  with  =  0  and 
e21  =  ®  1S  re9ar<ted  as  a  combined  expression  of  two  correlated  multivariate 
regression  models  with  different  design  matrices.  When  pj  =  p2  a  1,  Zellner 
(1962,  1963)  called  it  a  seemingly  unrelated  regression  (SUR)  equation  model, 
while  KFK  called  it  a  correlated  regression  equations  (CRE)  model.  As  has 
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been  discussed  above,  the  information  =  0  and  =  0  on  9  cannot  be 

expressed  as  a  single  GMANOVA  restriction  and  when  it  is  expressed  as  the 

» 

form  in  (1.5),  and  Mj  are  neither  nested  nor  orthogonal  relative  to  X  X 
(see  Kariya  (1985)  p  143,  159-163,  207).  Therefore,  the  model  itself  is  not 
put  in  the  framework  of  the  GMANOVA  model  and  even  the  LRT  for  Problem  til] 
is  difficult  to  derive  (in  the  case  that  XjX2  t  0  or  Xj  is  not  nested  to  X2). 

In  Section  .7,  based  on  the  idea  of  Mukherjee  and  Chandra  (1984),  we  compare 
the  power  of  the  LRT-like  test,  the  Pillai  type  test  and  the  LBI  test  proposed 
in  KFK.  The  comparison  is  made  asymptotically  in  n  with  contiguous  alternatives 
and  asymptotically  in  small  canonical  correlations  pj  ^...^pp  (as  p^O). 
There  it  is  shown  that  the  Pillai  type  test  is  asymptotically  (in  the  same 
sense)  equivalent  to  the  LBI  test  and  that  for  those  small  's  such  that 
t  =  lp.j/(ip..)  >  some  constant  c,  the  LRT-like  test  is  asymptotically  better 

than  the  LBI  test,  while  for  those  small  p^!s  such  that  t  <  c,  the  LBI  test  is 
better.  Further,  because  in  the  case  of  Xj  =X2  =  XQ,  the  model  in  (1.1)  is 
reduced  to 


[Y1,Y23  *  xo[eil*e22]  +  [E1'E2]» 

which  is  nothing  but  the  MANOVA  model.  Therefore  the  comparison  in  Section  7 
holds  as  it  stands  for  the  LRT  and  the  Pillai  (LBI)  test  of  independence  in  the 


MANOVA  model.  This  is  a  corollary  of  our  result. 
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2.  MOTIVATION  OF  THE  WORK 

The  motivation  behind  Problem  [l]:0j2  =  0  anc*  e2i  =  ®  1S  associated 
with  the  problem  of  no  causality  from  to  Y2  and  total  exogneity  of  Xj  for 
Yg.  We  will  first  write  the  model  (1.1)  as  two  correlated  classical  multi¬ 
variate  regression  models. 


Then  the  hypothesis  e12  =  0  is  equivalent  to  no  effect  of  Xj  on  Y£  as  in  the 
usual  case.  However,  since  Ej  and  E2  are  correlated,  the  regression  equation 
of  Y2  under  e12  =  0,  conditional  on  Yj  is  expressed  as 


y2  =  x2e22  +  (VXleirX2e2l)ZllZ12  +  e3< 


In  this  sense  the  effect  of  X^  on  Y-2  still  remains  unless  E^2  =  0.  Therefore 
the  hypothesis  e^2  =  0  and  lj2  =  0  in  Problem  [I]  is  considered  as  a  formu¬ 
lation  of  no  causality  from  Xj  to  Y2  or  total  exogeneity  of  X^  for  Y2>  An 
example  for  Problem  [I]  is  found  in  a  problem  of  economic  policy  evaluation. 
Suppose  the  model  (2.1)  is  a  reduced  form  of  an  econometric  simultaneous 
equations  model  which  describes  the  interaction  of  economic  variables,  and 
Xj  is  a  matrix  of  policy  variables  (tax  rate,  government  investment  etc.). 

Then  the  hypothesis  in  Problem  [I]  is  interpreted  as  no  effect  of  the  policy 
on  some  economic  variables  such  as  inflation  rate,  sales,  consumption  etc. 

The  motivation  behind  Problem  [II]  :.-6^ ^  =  02i  =  ®  anc*  Z\z  ~  ®  1S 

similarly  given  in  association  with  no  causality  from  Xj  to  Y2  and  no  causality 
from  X2  to  Yj.  In  addition,  the  problem  is  also  considered  as  a  formulation 
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of  no  additional  information  hypothesis  in  canonical  correlation  analysis, 
which  was  given  by  Fujikoshi  (1982)  based  on  Mckay  (1977).  To  see  this, 
suppose  there  are  two  groups  of  measurements  (variables),  say  and  x^,  where 
Xj's  are  r..xl  random  vectors  with  means  and  joint  covariance  matrix  y  = 

(^  j)  :(ri+r2)x(ri+r2)  with  4^ :  r^xr^.  (i,j  =  1,2).  Let  <52(xj,x2)  denote  the 
sum  of  squares  of  the  canonical  correlations  between  x^  and  x2 

52(xrx2)  =  tr^j  j4'12'lr22'i,2i  *  (2.2) 

which  is  regarded  as  a  measure  of  total  correlation  between  x^  and  x2>  Some¬ 
times  for  each  group,  there  are  some  other  measurements  available,  say  y ^  and 
y2,  which  appear  to  be  of  some  relevance  for  the  correlation  between  the  two 
groups  where  y^  :p^xl  ( i =1 ,2) .  Then  adding  these  variables  to  x^  and  x2,  the 
total  correlation  is  measured  by  the  sum  of  the  canonical  correlations  between 

21  =  ^  *  and  z2  "  (x2*yP'»  say  5 (zi*z2^  as  in  the  case  of  (2*2)*  But 

the  real  or  significant  relevance  of  including  the  additional  variables  y j  and 
y2  may  be  in  question  relative  to  the  original  variables.  This  question  gives 
the  following  testing  problem 

H:62(z1,z2)  =  62(x1,x2)  vs  K:62(z1,z2)  >  <S2(x1,x2).  (2.3) 

Using  a  conditional  argument,  Fujikoshi  (1982)  showed  that  this  problem  is 
equivalent  to  Problem  [II]  and  he  derived  the  LRT  where  X . ' s  and  Y.'s  in  (1.1) 

*  J 

are  the  sample  matrices  of  x. 's  and  y/s.  Mckay  (1977a)  treated  the 
2  2 

hypothesis  <5  (z^,x2)  =  6  (x^,x2)(i  .e. ,  no  additional  information  in  y1  relative 
to  (x^,x2))  and  showed  that  this  hypothesis  is  equivalent  to  in  the  model 
(1.1).  Some  related  topics  are  also  found  in  Mckay  (1977b)  and  Rao  (1970). 

The  motivation  for  Problem  [III]  is  stated  in  KFK  (1984),  Zellner  (1962, 
1963)  and  the  articles  therein. 


3.  TESTS  FOR  012  =  0  AND  Zu  =  0 


Based  on  the  motivation  stated  in  Section  2,  we  here  consider  the  problem 
of  testing  the  hypothesis  H  against  K  where 


H:  012  =  0,  =  0,  K:  not  H. 


(3.1) 


First  we  make  an  invariance  consideration  into  the  problem  and  obtain  an  expression 
for  the  local  behavior  of  the  power  function  of  an  invariant  test  in  the  neighbor¬ 
hood  of  the  null  hypothesis.  Based  on  the  expression,  an  invariant  test  together 
with  the  LRT  will  be  proposed  and  then  the  null  distributions  of  these  test  statistics 
will  be  given  in  Section  5.  To  begin  with,  a  canonical  reduction'crf“the  problem  is 


performed.  Write 


X  =  P[o]with  Pe°(n)  and  A  =  (X'X)‘  eG-e(r) 


and  express  0J2  =  0  as 


=  0  with  Mj  *  Hp^.O]  and  M2  =  (j 


(3.2) 


where  0( n)  denotes  the  group  of  n*n  orthogonal  matrices  and  Gf(r)  the  group 
of  rxr  nonsingular  matrices.  Further  let 


.-1 


MjA  =  F(I,0)f  with  FeG-^rj)  and  ¥  etf(r) 


and 


Z  = 


>  0 


0  I: 


P'Y  and  n  = 


n-r 


nn  ni2  \ri 


r  H'AQ. 


(3.3) 


n2i  ^22  | r2 


1  P2 


Then  the  problem  is  to  test  (1,0)  =  n19  =  0  and  E, 


=  0,  i .e. , 


in  the  canonical  model 


/zu 

Z12 

Z21 

Z22 

;Z31 

Z32 

Pi 

CM 

Q. 

r.  This  p 

nll  n12 

n21  n22 
0  0 


E11  E12 
E21  E22 


G  =  0(n)x8(p)xF  acting  on  Z  and  (n,r)  by 

g(Z)  =  PZB  +  F 


P1  0  \  F11  0  \ 

g(n,z)  «  (  I  nB  +  ,  B'eB) 


0  P, 


F21  F22 


(3.6) 

(3.7) 


where 


g  =  (P,B,F)eG 


P:  0  0 


0( n)  =  {P  =  0  P2  0  |  P.ecX^.)}, 

\o  pJ 


B(p)  =  {B  = 


B.e6£(pi)} 


F n  0  \.  r-jPi 

F  =  {  F  =  11  1|  F..eR  1  J} 

F21  F22 1 
\0  0  / 


It  follows  from  (3.7)  that  the  power  function  of  an  invariant  test  is  a  function 


of  (S12S12*n)  where 


^12  =  nl2^22  ^ ^  -  o(p)  =  jj  find 


(3.8) 


a  =  a(p)  =  diag{p1,...,pt}:  Pjxp2 


(3.9) 


2  2 

with  t  =  min(p1,p2).  Here  p1  ^  ...  >  pt  are  the  characteristic  roots  of 


-1-1  -i'  -1  -i' 

^ll^lZ^^l  an<*  $  1S  an  orthogonal  matrix  which  diagonalizes  I22L 21^11^12^22' 


Hence  without  loss  of  generality  we  assume  =  Sj2  ancl  Further,  it 

also  follows  from  (3.6)that  any  invariant  test  is  a  function  of  (Zj2,U)  with 
U  =  (Z31,Z32),  on  which  G  acts  by 


g(Z12,U)  =  (P1Z12B2,  P3UB)  for  g  =  (P,B,F)eG.  (3.10) 


Now  to  state  one  of  our  main  results  in  this  section,  let  VQ  be  the  set  of 


all  invariant  tests  of  size.*. (Ke  v  $ efr- <->  fC?(Z))  =  <j> ( Z) ) , 


S  -  U'U  5  L11  •  Sij  =  Z3iZ3j 


S21  S22 


(3.11) 


and  let  6  = 


+  52  with 


,  '  -l  ' 

61  =  tr512512  =  trr,12E22n12 


S2  =  l  i=lpi  =  trEllE12Z22Z21 1 


(3.12) 


Clearly  n^2  =  0  and  e^2  =  0  if  and  only  if  6=  0,  or  6^  =  0  and  S2  =  0. 


Theorem  3.1  .  There  is  e  >  0  such  that  on  the  set  { (n ,E ) 1 6  <  e),  the 


power  function  of  any  test  <j>  in  V *  is  evaluated  as 


tt  (4>  ,(n,z))  =  a  +  <5^Cj(<j>)  +  2C  2  ( )  +  0(5) 


(3.13) 


C1(<()) 


^  <  rj  E0^trZ12Zr/Z12Z12+S22^  »•  <3*14> 


c2U) 


rl+r3  c  r 


.-1 


2pxp2  E0{‘i,tr3trSllS12^S22+Z12Z12)  ^21 


(3.15) 


P1trS22(S22+Z12Z12)  }  -  23, 


lim  supx  I o(6)/6 I  =  0. 

6+0  + 

Here  Eq(0  denotes  the  expectation  of(*)imder  the  null  hypothesis.  The 
proof  is  given  at  the  end  of  this  section.  The  expression  (3.13)  shows  the 
local  behavior  of  the  power  function  of  an  invariant  test  <|>,  according  to 
which  the  power  function  is  approximated  by  a +  6^(41)  +  62C2(<|>)  in  the 
neighborhood 

-1  1  -1  -1 

N£  —  ( (n  ,E )  ]  trn ^2^22^12  ^^*^11^12^22^21  ^ 

and  the  approximation  is  uniform  in  <j>.  Since  a  test  maximizing  the  local  power 
a  +  6^Cj(<j> )  +  §2  C2 (<t> )  depends  on  at  least  the  ratio  of  6j  and  62  ,  no  LBI 
test  exists.  This  is  natural  in  the  sense  that  the  quantities  6j  and  62  indicate 
are  different  and  the  deviation  of  6j(or  n^2)  from  0  is  independent  of  the  devia¬ 
tion  of  62  (or  Zj2)  from  0.  On  the  contrary,  the  form  of 

the  power  function  in  £3.13)  reflects  how.  an  invariant  test  can  detect  each  local 
deviation  from  each  null  case.  That  is,  Cj($)  basically  measures  the  local  power 
of  <p  against  the  local  deviation  of  6j  from  0,  while  C2U)  the  local  power  of  <(> 
against  the  deviation  of  52  from  0. 

However,  the  statistics  whith  define  the  expectations  of  0.(41)  and  C_(4> ) 


in  (3.14)  and  (3.15),  are  not  independent  even  under  the  null  hypothesis.  In 
fact,  the  statistics 


s  rj  trZ12Z12^Z12Z12+S22^  1  and 

-  r3trSllS12^S22+Z12Z12^  ls2l"PltrS22(S22+Z12Z12) 


(3.16) 


are  dependent  on  each  other.  This  is  a  feature  of  the  simultaneous  treatment 
of  the  two  seperate  hypotheses.  To  investigate  this  point  further,  observe 
that  the  test,  say  ipj,  which  maximizes  Cj(<j>)  is  given  by  the  critical  region 

T^  >  c 

and  it  is  the  LBI  test  for  testing  =  0  without  =  0*  More  specifically 

it  is  LBI  for  testing  the  General  MANOVA  hypothesis 

H ' :  Mj  ©Mg  =  0  (3.17) 

in  the  MANOVA  model  Y  =  XQ  +  E  with  E  ~N(0,I  ®  e)  where  X:  n*r,  M^:  r^xr  and 
M2:  pxp2  are  arbitrarily  fixed  matrices  of  full  rank.  This  test  is  even  UMPI 
(uniformly  most  powerful  invariant)  when  r^  =  1  or  p2  =  1  and  the  power  function 
of  an  invariant  test  *l>  for  the  hypothesis  in  (3.17)  is  locally  expressed  as 

TT1(i'  »(n,z))  =  a  +  SjCjOP)  +  ©(Sj)  (3.18) 


where  (n,E)  is  the  parameter  of  a  canonical  form  corresponding  to  (e,E)(see 
Kariya(1985)  ,  p  109).  On  the  other  hand,  the  test,  say  t|;2,  which  mazimizes 
C2 ( 4> )  is  given  by  the  critical  region 


with  in  (3.16)  and  it  is  the  LBI  test  for  testing  e^  =  0  in  the  case  of 

n12  =  ^ore  specifically  it  is  LBI  for  testing  independence  E^  =  0  in 
the  missing  data  model 

Zj2  ^^12*^r  ^22^  with  n^2  —  ^ 

1  (3.19) 

U  '  N(0,Ir  ®  z)  and  and  U  are  independent 
3 

where  the  counterpart  of  ZJ2  is  missing  (see  Eaton  and  Kariya  (1983))  and 
the  power  function  of  any  invariant  test  in  the  model  (3.19)  is  expressed  as 

ir2(^»(0,E))  =  a  +  62C2^  ^  +  °^S2^  (3-2°) 

However,  when  n12  t  °»  the  test  based  on  the  critical  region  T2  >  c  is  not 
easy  to  interpret  as  a  test  for  testing  the  single  hypothesis  of  the 
independence  e12  =  0  alone  because  in  the  case  of  n12  *  °»  Z12  would  not  be 
involved  in  a  test  statistic.  In  other  words,  this  is  a  difference  between 
treating  the  two  hypotheses  simultaneously  and  treating  them  separately. 

Now  by  taking  this  point  into  account,  we  propose  the  test  maximizing 
a  linear  combination  of  Cj(<J>)  and  Cg^); 

C  U)  =  ^r3C1  (4> )  +  C2U)  (0  <  6  <  »)  (3.21) 

where  B  is  a  constant  independent  of  n.  Using  the  Generalized  Neyman-Pearson 
Lemma,  the  critical  region  is  given  by 

T(B)  :  Sr^  +  T2  >  k.  (3.22) 

Here  is  multiplied  by  r^  =  n-r  because  from  (3.3)  n  =  0(n*)  so  that 
61  =  0(n)  provided  X'X  =  0(n).  The  test  <j>g  with  critical  region  T(e)  >  c 
maximizes  the  power  irU,(n,z))  locally  in  the  neighborhood 


since  It  maximizes  (Sj/r^Br^U)  +  62C2($)  =  figleir^U)  +  C2 (<i> ) ]  on  Neg. 

The  constant  6  may  be  regarded  as  a  weight  for  the  importance  of  the  hypothesis 
n12  =  0  re1ative  to  the  hypothesis  s12  =  0,  and  it  is  chosen  in  advance.  It  is 
noted  that  the  test  based  on  T(b)  is  not  a  linear  combination  of  the  two  LBI  tests 
and  ^2  stated  above. 

We  remark  that  for  a  given  <j>  e  P*,  the  local  sensitivity  of  $  against  (6j,62) 
is  measured  by  the  two  coordinate  (C^(4»)  ,C2(<^>)) .  Second ,  when  the  information  ri^2  =  0 
is  ignored  in  the  missing  data  model  (3.22),  the  LBI  test  for  independence  E12  "  0 
is  given  by  the  critical  region 

T3  5  trSjjsl2S22S21  >  k  (3.23) 

(see  Schwartz  (1967)).  That  is,  this  is  the  LBI  test  for  independence  in  the 
MANOVA  model  without  n12  =  0.  Since  the  test  with  critical  region  Tj  >  k  is  LBI 
for  testing  n12  =  0  without  lj2  =  0,  we  may  combine  these  two  tests.  Here  while 
the  simultaneous  treatment  of  the  hypotheses  yielded  the  dependent  statistics  Tj 
and  T2,  the  separate  treatment  yields  the  independent  test  statistics  Tj  and  T3. 

The  problem  of  how  to  combine  independent  test  statistics  is  discussed  in  the 
literature  (see,  e.q.,  Marden  (1982)).  In  this  paper  we  do  not  get  involved 
in  this  problem.  But  as  is  shown  next,  the  LRT  for  the  simultaneous  hypotheses 
gives  a  natural  combination  of  the  LRT  statistics  for  two  separate  hypotheses. 

The  LRT  for  our  problem  is  easily  obtained  by  using  the  canonical  model  in 


(3.5)  as  follows: 


Of  course,  when  is  small,  the  joint  hypothesis  is  rejected.  As  is  well 
known,  Lj  is  the  LRT  statistic  for  =  0  without  =  0  while  is  the  LRT 
statistic  for  r^2  =  0  without  lj2  =  0.  Further  Lj  and  L2  are  independent  so  that 
the  two  independent  LRT  statistics  are  combined  in  (3.24).  It  is  noted  that 
even  when  n^2  =  0»  the  LRT  statistic  for  Zj2  =  0  is  given  by  L^  so  that  the  LRT 
ignores  the  additional  data  (see  Eaton  and  Kariya  (1984)).  The  unbiasedness 
property  of  the  power  function  is  considered  as  a  special  case  of  Problem  (II) 
in  the  next  section. 

Proof  of  Theorem  3.1.  The  proof  is  similar  to  KFK  (1984).  To  derive  the 
distribution  of  a  maximal  invariant  T  =  T  =  T(Zj2,U)  under  the  action  (3.10)  of 
G  on  (Z12,U),  we  apply  the  Wijsman's  representation  theorem.  Let  pj^  ^  be 
the  distribution  of  a  maximal  invariant  T.  Then  the  density  of  T  with  respect 
to  P^0  evaluated  at  T  *  T(Zj2,U)  is  given  by 

R  -  dPU12,n)/dP(0,I)  s  H(Z12,U|c12,o)/H(Z12,U|0,I)  (  3.2$) 

where 

H( zi2 ,11  !5i2 ,«)  -  f  f(PiZ12B22’  P3UBU12,n)X(B)v1(dP1)v3(dP3)y1(dB1)y2(dB2) 

•  H 

■  *3/2  1  (^*i  )/2  i  -Pi/o 

X(B)  =  |B1B1|J  |B2B2|  1  J  ,  Uj(dBj)  =  IBjBjI  J^dBj 

H  =  0( rj)  x  G(r3)  x  {^(p^  x  &d(p2),  f (z12 ,U | C12 »«)  is  the  density  of  ZJ2  and 
U  with  Z12  ~  N(c12,I®  I)  and  U  ~  N(0,Ir  ®n),  and  v^dP.)  is  the  invariant 
probability  measure  on  0(r.)  ( i =1 ,3 :  j=l,2).  The  condition  for  which (3.25 ) 
holds  is  satisfied  (see  Wijsman  (1967)).  In  order  to  obtain  the  local  behavior 
of  the  power  function  of  an  invariant  test  in  a  neighborhood  of  C^2  =  0  and 
z,«  =  0,  we  evaluate  R  in  (3.25  )locally.  After  cancellation  of  some  constants  the 


numerator  in  (3.25)  is  expressed  as 


H(Z12.U|c12,fi)  =  j 


_  p 

M  3/2  exp[-J5K]X(B)v1(dP1)y1(dB1)u2(dB2) 


where 


K  -  ^r^2^12^12®2  "  ^^1^12^2  h.2  +  tr^12^12 


(3.26) 


I  I 


+  trBis1iB1F11  +  2trB1S12B2F12  +  trB2S22B2F22 


with 


.-1 


I  c  r 

rll  h12 
F21  F22 


(3.27) 


Here  replacing  Bj  by  sih  and  B2  by  (zj2^12+S22^”35B2  leaves  the  ratl"°  R 
remain  the  same  since  the  Jacobin  factors  coming  out  are  cancelled  out  with 
those  of  the  denominator.  Further  writing  F^  as 

Fii  =  tI‘A1i)  1  =  1  +  An  +  Aii(I'Aii}  1 

I  I 

with  =  AA  and  a22  *  a  a.  Then  K  in  (3.26)  becomes 
K  =  trBjBj  +  trB2B2  -  ZtrPjZ^B^^  +  tr512q2  +  2trBlS12F12  +  Ki  +  K2 


K,  »  trBjBjlAj+Ajd-Aj)’1];  K,  =  trB^S„B9[A„+Ai(I-A,r1] 


2  22  2 1  2  2V  2J 


Z12  =  Z12^Z12Z12+S22^  *’  ^12  =  S11S12^Z12Z12+S22^  *  and 


where 


.  ’<-]  "•r,1  j  rj 


S22  =  (Z12Z12+  S22)  *S22^Z12Z12+S22^ 


Now  expand  the  components  of  exp[-y<]  as 


exp[trPjZj2B2c12]=  1  +  tf*P 2^X2^2^12  +  ^^|,'^1Z12^2^12^ 


•  -  i 


■  -  i 


exp[-trB1S12B2F12]  =  1  -  trB1S19B9F19  +  %(trB1S1  ,B,F., +o 


1122  12 


'ri2u2'12; 


expt-JsKj]  =  1  -  J^trBjBjAj  +.  Og  and 
expt-J^]  =  1  -  %trB2S22B2A2  +  c>4. 

Then  in  the  same  way  as  in  KFK(1984),  the  remainder  terms  o^'s  are  shown  to  be 
e(s)  with  <5  in  (3.12),  which  is  uniform  in  (Zj2,S)  because  of  the  boundedness 
of  Z12  and  S^'s  .  Further  when  o^'s  are  integrated  over  Pj  and  B^ 's,  they  are 
shown  to  be  o(6)  uniformly  in  (Z12,S).  Now  taking  the  product  of  these  terms 
and  ignoring  the  odd  functions  of  eitherPj  or  Bj  or  B2  because  they  are  zero 
when  integrated,  we  obtain 


R  =  |fi|  r3/2J  Jv1(dP1)h1(B1)h2(B2)dB1dB, 


where 


h^Bj)  =  exp(-%trB.B! )  |B.b!  |^i/2dB./D. 

J  =  1  +  ( 'tv'P ^ “  trB1B1  AA 


1122  12J 


t  -  I 


-  trB2S22B2A  A]  +  Og  =  1  +  ^[I+II+III+IV]  +  o5,  say 
0i  =  exp(-J5trBiB!)|B.Bj!Mi/2dBi, 


.‘•V  .  -  V-  .>V-  ‘  *  r"  -  ■  -  \  \-  *.*  v  V  *  'V  v.-  .*  .-\  V  V* 

*  %“•  %*.  V  */  *.\V  .  .  w*..  „ a  *-  V  ,*•  --  ^ 


M1  =  r3  "  P1  and  M2  =  rl  +  r3  "  p2*  Here  H(Zk,U(0,I))  =  and  |n|  3^2  = 

1  +  o(<s)  were  used.  The  remainder  o,.  is  shown  to  be  o(6)  for  6  small  uniformly 
in  (Zj2»S)»  for  =  (I-AA  )  A.  We  now  need  to  integrate  J.  But  arguing  as 
in  KFK  (19S4)  pp  391-392,  we  obtain 


I  = 


II 

III  = 

IV  = 


{trP1Z12B2?12)2v1(dP1)h2(B2)dB2  -  tr^12C12tr^12^12 ! 

-  |(trBjs12B2F12)2h1(B1)h2(B2)dB1dB2  =  tr£j2S12trAA  +  0(6); 


1  2 

ii  ^3^1  1 

(trB^BjAA  Jh^B^dB^  = - p —  trAA  , 

1 

(trB2S22B2A  A)h2(B2)dB2  = - p — 3  trS^trAA 


and  |0g  j(dPj)hj(Bj)h2(B2)dB^dB2  =  o(&) .  Thus,  observing  that  the  power  function 
of  an  invariant  test  <j»  is  given  by 


*u.(n.r))  =  J+dP(e12n)  *  /*Rdp(o.i)* 


4.  TESTS  FOR  THE  HYPOTHESIS  012=O,  6^=0,  AND  Z12=0 


We  shall  first  consider  via  invariance  the  Problem  (II): 


H:  0j2  =  0,  =  0,  Xj2  =  0  vs  K:  not  H.  (4.1) 

Since  e12  =  0  and  e21  =  0  are  expressed  as 


M^eM^  =  0  and  M^gM^  =  0 

I  I 

and  since  M^X  XM3  =  X^X2  f  0  in  general,  the  group  leaving  this  problem  invariant 
is  smaller  than  the  group  leaving  the  problem  (I)  invariant,  where  and  M2  are 
defined  by  (3.2  ),  Mg  =  (0,1)  and  M^  *  ^  .  The  special  case  X^X2  =  0  will  be 
briefly  treated  later.  Here  we  use  the  following  canonical  form 


W  =  K(x'x)‘1x’y  *  N(n,A<g>z)  with  o  -  Kg 

S  *  y'(I-P0)Y  “  W(z,r3)  with  r3  =  n  -  r 
W  and  S  are  independent 


where  Pn  = 


x(x'x)"1x'. 


K  = 


I  q;^  0  ' 


<11 
0  ,  Q 


-*s 

22 


with  (x  x)-1  -  q  -  f Ql1  °12ri 

tQ2i  $22  1' 2 


rl  r2 


and 


A  =1 


I  0 

iq;2  i 
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with  Qj2  =  Qj^Qj2^22  ’ 


(4.2) 


(4.3) 


(4.4) 


Partition  W  and  n  as 


■I' 


J"  11  "12  '  1 

\ W21  W22  r2 
'  n  n  t 


and  n  = 


"11  "12  '1 


n21  n21  /  2 
Pi  P2 


respectively. 


P1  p2  P1  p2 
Then  in  the  model  (4.2)  the  problem  is  to  test 


H:  n12  =  0,  n21  «  0,  zu  =  0. 


rlPl  nr2P2 


The  problem  is  clearly  left  invariant  under  the  group  G  =  G£(Pi)*G£(P2)*R  *R 


acting  on  (W,S)  by 


,  B,  0  F,  0 

g(W,S)  =  (WB+F.B  SB)  with  B  =  *  R  and  F  =  nA  . 


(4.5) 


where  g  =  (Bj.B^F^F^eG.  From  this  action  of  G,  it  is  easy  to  see  that  an 
invariant  test  is  a  function  of  (W12>W21>S)  only  and  that  the  power  function  of 
an  invariant  test  is  a  function  of 

n  s  J,  A  '  512  “  Q11  012I:22  P2  and  KZl  *  QZ2QZlZn  P1  (4'6) 

r  l\ 

where  P.etfp^'s  (i=l,2)  satisfy 

PlEilEl 2Z22P2  =  A  =  diag^  pi»*  pixp2  with  1  =  min(Pj  *p2^ 

and  p  >  >  pf  are  the  canonical  correlations.  So  we  can  assume  n12  =  Cj2» 

1 

n  =  5^  and  E  =  n  without  loss  of  generality.  In  this  set-up,  we  may  proceed 
in  the  same  way  as  we  did  in  Section  3  to  obtain  an  expression  for  the  locally 
approximate  power  of  an  invariant  test.  However  by  doing  so  we  end  up  with  a 
very  complicated  expression  of  local  power  which  depends  on  many  parameters  in  an 
intractrable  way.  This  implies  that  a  further  invariance  consideration  into  the 
problem  will  not  help  to  propose  a  new  test  which  possibly  takes  into  account  the 


22 


simultaneous  occurance  of  the  three  hypotheses  n12  =  0,  n21  =  0,  and  z  =  0.  Of 
course  no  LBI  test  exists. 

On  the  other  hand,  the  LRT  is  given  by  the  critical  region  L  <  c,  which 
was  originally  derived  by  Fujikoshi  (1982),  where 


^1^2^3*  sa^ 
(4.8) 


It  is  also  directly  obtained  from  the  distribution  of  (Wj2,W21,S)  since  the 
LRT  is  always  invariant  (under  a  very  mild  condition).  The  statistics  Lj,  L2 
and  in  (4.8)  are  respectively  the  LRT  statistics  for  the  three  separate 

hypotheses  =  °>  e12  =  and  021  =  °*  and  the  LRT  statistic  LiL2L3  for  our 
simultaneous  hypothesis  in  (4.8)  may  be  viewed  as  showing  how  to  combine  the 
three  LRT  statistics  for  the  three  separate  hypotheses.  But  here  it  is  noted 
that  L2  and  L3  are  correlated  because  of  the  correlation  of  W12  and  W21  unless 

I  II 

XjX2  =  0  or  E12  =  0.  In  fact,  if  u.  and  are  respectively  the  i-th  and 
j-th  rows  of  Wj2  and  W21,  the  covariance  matrix  of  un-  and  v^.  is  shown  to  be 

I 

Cov  (i^.Vj)  =  q^.A  :  p^Pj 


where  q^.  is  the(i,j)th  element  of  $12  in  (4.4)  and  fJ12  -  0  and  A  =  0  are  res  pec- 

I 

tively  equivalent  to  X^X2  =  0  and  =  0.  This  correlation  between  L2  and  L3 
makes  it  d.rficult  to  investigate  optimality  properties  of  the  LRT.  Now  to  show 
the  unbiasedness  of  the  LRT,  note  that  from  wj2'NUj2»I  ®  I)  and  W2l'N^21’*®  I) 

I  * 


^22  *  ^12^12  (I jT2:t2^  with  ^2  “  Pj  and 


=  C12512 


(4.9) 


=  w2iw2f  wpi(I»ri:Ti) 


with  p^  “ r 2  and  *  ^2X^21 


-  '  *'•**'  ■.  *  L  ■  - 


b.  V..  *" 


w  »  .  " 


>  V 


."vVv 


>.  V 'I- 

■  .--1  jJ 


where  W('i',m:v)  denotes  the  noncentral  Wishart  distribution  with  mean  m^,  degrees 


of  freedom  m  and  noncentrality  parameter  \>.  In  the  original  term,  V..  i 


shown  to  be  equal  to 


vii  =  VVpi]Yi  with  pi  =  xn- <xixi >"lxi  -  (4-10) 


Further,  from  (4.2), S  and  {V^,V221  are  independent,  and  when  A  =  0, 
are  independent.  Here  we  use  a  conditional  argument.  First  write 
(4.8)  as 


L1L2  =  lS22-ll/lS22-l  +  V22  +  S21S11S12  ^  =  ^22*1^/^22*1  +  ^22' 


where 


S22 - 1  =  S22  ‘  S21S1 1S 12 *  ^22-1  =Y  ^  S22 -1T 


U22  *  y’,*(V,9+S,1s;}s1,)y’,s  and  y  *  I  -  a’a. 


22  21  11  12 


Lemma-  4.1  (1)  S22l~  W  (I.n-r-pj) 


(2)  S22  ^nd  U22  are  independent 


(3)  Conditional  on  Yj,  U22~  Wp  (I ,r1+P1 :f)  where 


'v  »  Y  '^[e+a'S^Iy-*5 


(4.12) 


[^1^12+^1^]  ^0*^2^  with  ^2  ”  ^12^22^2 


Proof.  Our  original  model  may  be  viewed  as 

0  6 . 


[YrY2]  '  N((X1,X2] 
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Le21  0  ^ 


,  in  ®  n) 


with  n  in  (4.6),  where  e?,  =  0?1E_is-P, .  Hence  conditional  on  Y. , 


easily 


and  V22 
in 


ID 


Y2  '  N(X1012+{Y1-X2821)A,  I  ®  Y) 


from  which  conditional  on  Y^,  we  obtain 

S21S11S12  "  Y2P0V2  '  Wp2^’pl:4  Sll4'’ 

V22  =  V>0-P1!V1  '  “p2(’'-rl:5)  an<l 
S22-l  ’  “p2W.n-r-Pi). 

Further  conditional  on  Y^,  ^22  •  1  ’  $21^11^12  anc*  Y22  are  indePendent  and 
S22  ^  does  not  depend  on  Yj.  Thus  all  the  result  follow. 

Theorem  4.1.  The  LRT  with  critical  region  L  <  c  is  unbiased. 

Proof.  Let 

C={(S22*1’  V22+S21S11S12’V11,S11):  L1L2L3  <  C} 

be  the  critical  region  of  size  a  in  the  space  of  (s22-l’  V22+S21SllS12’VirSll^ ' 
Then  since  Sjj  and  are  functions  of  Yj  only,  from  (4.10)  the  power  function 
of  the  LRT  is  expressed  as 

if ( 4* » (® ^2 *® 21  ^  -  I ®  12  *®21  ^ 

=  Ey^  [P(  { (^22-i  »^22)  I  1^22*1^^22*1  +  ^22^  <  C^3  ’^11  *®12  *®21 (4.13) 

where  Ev  denotes  the  expectation  with  respect  to  Y., .  Since 


S22*l  +  U22 I 


>: 

v. 


25 


is  regarded  as  the  LRT  for  testing  V  =  0  in  the  MANOVA  set-up  with  (U22,S221), 
as  is  shown  in  Anderson  et.al.  (1964),  it  is  an  increasing  function  of  each 
characteristic  root  of  f.  Since  §12  =  0  and  A  =  0  imply  f  =  0,  from  (4.13) 

^  21  ))  *  P(  C  |  0 ,0  2^  1 0 ) .  (4.14) 

But  under  A  =  0,  L^L2  and  are  independent  because  and  V22  are  independent. 
Further  =  IS^I/IS^  +  is  regarded  as  the  LRT  statistic  for  testing 

521  =  0  1n  the  MANOVA  set-up  with  (V^.S^).  Hence  the  inside  of  the  conditional 
expectation 

p(c|o,§21,o)  =  el  l  ip({(s11.v11)|l3  <  c(L1L2)'1}|(O,021,O)] 

w  I 

is  an  increasing  function  of  each  characteristic  root  of  ©21® 21  ^mP"' ”*es 

P(CjO,021,O)  >  P(C|0,0,0)  >  a 

Combining  this  with  (4.13)  yields  the  result. 

Under  a  t  0  or  equivalently  i^2  f  0,  because  LjL2  and  are  correlated 

I 

unless  XjX2  =  0,  it  seems  difficult  not  only  to  establish  a  monotonicity 
property  of  the  power  function  of  the  LRT  but  also  to  find  a  naturel  parameter 
space  on  which  the  monotonicity  is  considered. 

We  remark  that  the  above  result  holds  even  when  the  model  Y  =  Xe  +  E  is 
defined  for  X  conditioned  and  the  marginal  distribution  of  X  does  not  depend  on 
(e,z) .  Hence  the  LRT  for  testing  no  additional  information  hypothesis  in 
canonical  correlation  model,  which  is  nothing  but  our  LRT  though  X  is  random, 
is  unbiased  (see  Section  2). 

I 

In  the  case  of  X^X2  +  0,  we  consider  two  special  cases.  First  consider 


the  case  A  =  0  or  equivalenty  e^2  *  0.  In  this  case  the  problem  of  testing 
joint  hypothesis  812  =  0  and  e21  =  0  is  simply  split  into  two  independent 
problems  in  the  two  independent  models  Y..  =  X^e +  X202l-  +  E^(i=l,2).  However, 
the  LRT  for  the  joint  hypothesis  is  given  by  <f>L  =  X( L2L3  <  c)  and  it  is  not  the 
product  of  the  two  LRT  X(L2  <  c2)X(L3  <  c^)  where  X(A)  denotes  the  indicator 
function  of  a  set  A.  The  power  function  of  <J> ^  for  0^2  =  0  and  0^  =  0  is  a 
function  of  the  characteristic  roots. 

=  ch,.  ^22®  12^1  ^0~^2  12^  22  ^  •  •  »P£) 

and  the  characteristic  roots 

Yj  ~  ch j ^11L^ 21  ^2  ^0”^1  ^2^ 21^1 1  ^  ( j=l » •  •  •  »Pj ) » 

and  by  similar  argument  as  in  the  proof  of  Theorem  4.1,  we  obtain 

Corollary  4.1.  Under  E^2  ■  0,  the  power  function  of  the  LRT  with  critical 
region  L2L3  <  c  is  increasing  in  each  X.  or  Yj* 

Next  we  consider  the  unbiasedness  of  the  LRT  with  critical  region  L^L2  <  c. 
Since  the  distribution  property  of  LjL2  is  simply  obtained  in  the  proof  of 
Theorem  4.1  by  setting  V-q  =  0,  we  obtain 

Corollary  4.2.  The  LRT  with  critical  region  L^L2  <  c  for  Problem  [ I ] 
is  unbiased. 

In  fact,  from  (4.10)  and  Lemma  4.1,  the  power  function  of  the  LRT  is  easily 
seen  to  be  an  increasing  function  of  each  characteristic  root  of  f  conditional 
on  Y.  Hence  the  unbiasedness  immediately  follows  from  y  ^  0  and  the  fact  that 
e12  =  0  and  E12  =  0  imply  Y  =  0. 

Finally  we  consider  Problem  [II]  under  X,X„  =  0.  An  example  for  this  case 


is  found  in  time  series  models  (see  Anderson  (1971  p. 92)).  Under  =  0,  we 

can  take  A  =  I  in  the  canonical  form  in  (4.2)  as  well  as  Q..  =(X.X.)-1.  Then 

n  i  r 

the  group  G  =  0(r)xG(r2)*G  leaves  the  problem  invariant  by  the  action 
g(W.S)  =  (tWB+F.b'sb)  with  r  = 

where  g  =  (P^.Pg.g)^.  Hence  in  the  same  way  as  above,  the  problem  is  reduced 
to  the  problem  of  testing  =  ^1  =  0  and  A  =  0  in  the  canonical  form 

Wj2  N(Ci2»  I  ®  1)  with  =  ^11^12^22^2 

<  W21  '  N(c21,  I®  I)  with  C21  =  Q22621Ziiri 

S  “  W(n,r~)  with  n  in  (4.6  ) 

V  J 

where  Wj2,  W21  and  S  here  are  independent.  Also  G  acts  on  (V^.W^.S)  by 

g(W12,W21,S)  =  (PjW^Bg.PgWgiBj.B'sB). 

Since  the  group  G  is  bigger  than  G,  a  result  corresponding  to  Theorem  4.1  can 
be  derived.  To  see  this,  let  v  =  +  v2  +  with 

V1  =  trC12C12  =  trXlXl012Z22012 

'  1  -1  1 
Vg  ~  ^r^21^21  =  ^^^2^2^21^11^21 

v3  =  ^i=lpi  =  trAA  =  trIllE12I:22r21 

Theorem  4.2.  There  is  an  e  >  0  such  that  on  the  set  {(g,e|v  <  e)  the 
power  function  of  an  invariant  test  4>  of  size  a  under  G  is  evaluated  as 


where 


DjU)  =  Eq  [4>R^  1 


rN+r,  P-i  i  i  _i 

R1  =  FJ>y  rj'  trW12W12(W12W12+S22^ 


r j+r ~  p«  i  i  _i 

R2  =  TJP^  r^  trW21W21(W21W21+Sll) 


R3 


( rl+r3 ) ( r2+r3 ) 


>-l. 


P1P2 


tr(W21w2i+sn)  si2(W12W12+S22)  S21 


rl+r3 


k-l  r2+r3 


Pj  trS22^W12W12+S22^  ‘  P2  trSll^W21W21+Sll^ 


-1 


and  lim  sup  |o(v)/v|  =  0. 
v->0  * 

Proof.  The  proof  is  completely  similar  to  that  of  Theorem  4.1.  Regarding 

Z12  aS  W12*  rePlace  Bi  by  W21^W21W21+Sll^_l5Bl ’  Then  every  step  9oes  throu9h 
and  the  result  is  obtained. 

In  this  expression,  a  symmetry  which  is  lacking  in  Theorem  4.1  is  secured, 
and  the  statistics  Rj,  R2  and  R3  are  respectively  the  LBI  tests  for  testing 
the  separate  hypotheses  H^:?^  =  0,  H2:521  =  0  and  H3:i12  =  0  under  £12  =  0  and 
?2i  =  0  (see  Eaton  and  Kariya  (1983)  for  H3)  Also  for  each  hypothesis  ,  the 
power  function  of  an  invariant  test  is  expressed  as  o  +  +  o(v.)(1=l,2,3). 

Here  again  it  is  noted  that  the  statistics  Rj.Rg  and  R3  are  not  independent 
under  the  null  hypothesis.  Rj  and  R2  are  independent  under  the  null  hypothesis. 
Here  following  the  discussion  in  Section  3,  we  may  propose  a  combined  test  of  these 
statistics  with  critical  region 


s  s 


m 
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RCBj.Bg)  *  8ir3Ri  +  82r3R2  +  R3  >  c  (4.15) 

where  r3  is  put  on  R^'s  because  \>^  ■  0(n)  (1-1,2).  Here  the  constants  8^'s 
may  be  considered  indicating  relative  weights  for  the  three  hypotheses  ■  0* 
e21  *  0  and  z12  *  0-  It  will  be  often  the  case  *  e2.  This  test  clearly 
maximizes  the  local  power  a  +  In  the  neighborhood 

{(0,S)|  v.  =  r381.v3(1»l,2)}n  { (0  ,E )  |  v  <  e> 

Also  the  local  sensitivity  of  an  invariant  test  $  against  (vj,v2,v3)  Is  described 

by  the  coordinates  (Dj($) ,D2(<fr) ,  03($)).  Further,  from  the  observations, 
above,  we  may  use  the  test  based  on  R( B2)  in  (4.15)  even  if  xjx2  +  0. 

It  is  remarked  that  the  LRT  statistic  in  this  case  is  of  course  the  same 

I 

as  L  s  LjL2L3.  But  here  because  XjX2  *  0,  in  addition  to  the  independence  of 
Lj  and  L1 (1*2,3),  L2  and  l3  are  independent  though  L^,  L2  and  L3  are  jointly 
dependent . 


5.  ASYMPTOTIC  NULL  DISTRIBUTIONS  OF  THE  TEST  STATISTICS 


The  asymptotic  null  distribution  of  the  LRT  for  Problem  [II]:  e12  =  0, 
e2j  =  0  and  =  0  has  been  derived  by  Fujikoshi  (1982)  in  the  context  of 
the  problem  of  testing  no  additional  information  hypothesis  in  canonical 
correlation  anaylsis.  From  a  more  general  viewpoint,  we  here  briefly  treat  it 
in  a  systematic  way  and  then  consider  the  asymptotic  null  distribution  of  the 
LRT  for  Problem  [I]:e12  =  0  and  i12  =  0.  The  notation 

*  ~  Ap(q,n) 

denotes  that  for  independent  Wishart  matrices  A  and  B, 


X  =  | A | / | A  +  B |  with  A  -  W  (E.n)  B  ~  W  (l,q) 

K  K 

Lemma  5.1  Let  X^  ~  Ap  (q^  .n-d^)  and  X^'s  be  Independent  (i=l ,2).  Then 
P(-m  log  XjX2  <  x)  =  Gf(x)  +  [Gf+4(x)  -  Gf (x) ]  +  0(m"3),  (5.1 


where  G^(x)  is  the  cdf  of  x  vf) 


fl  *  plql’  f2  =  p2q2’  f  =  fl  +  f2’  m  =  n’p 

p  =  7  {fi[dr  \  <VP1'1)]  +  f2td2‘  l(vp2_1)]} 

v  s  s  +  [^i-d2+  -j(P2-p2-qi+cl2^^ 

s  =  Sj  +  s2,  Sj  =  -ft  (pj+q^-5)  and  s2  =  (p2+q2-5). 


Proof.  The  result  follows  directly  by  the  usual  method  based  on 
characteristic  function. 

Now  for  Problem  [II],  the  LRT  statistic  is  given  by  L^L^  in  (4.8)  and 
from  Lemma  4.1  and  (4.11),  it  is  easy  to  see  that  under  the  null  hypothesis 
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X2  5  L1L2  '  Ap2^rl+Pl*  n_r'Pl^  (5.3) 

X1  =  L3  ~Ap1^r2+P2*  n_r^  (5.4) 

and  x^  and  x 2  are  independent.  Therefore  the  following  result  follows  from 
Lemma  5.1. 

Proposition  5.1.  For  Problem  [II] :  0j2  =  0,  Q^\  ~  0  and  ~  t*ie 
asymptotic  null  distribution  of  LRT  based  on  -m  log(LjL2)Lg  is  given  by  (5.1) 

where  in  (5.2)  q^  =  r2  +  p^,  ^  =  rl  +  P1 ’  dl  =  r  and  d2  =  r  +  p. 

On  the  other  hand,  for  Problem  [I],  the  LRT  statistic  is  given  by  LjLg 
in  (3.24),  which  is  the  same  as  the  LjL2  in  Problem  [II]  and  hence  from  (5.3) 
under  the  null  hypothesis 

x2  *  L1L2  *  Ap2(rl+Pl-n-r-':,l’ 

The  asymptotic  null  distribution  of  ^2m  log  X2  is  well  known  . 

Next  we  consider  the  asymptotic  null  distribution  of  R(B^,B2)  in  (4.15) 
for  Problem  [II]  and  as  a  special  case,  obtain  the  null  distribution  of  T(b) 
in  (3.22).  But  here  for  simplicity,  the  case  B^r.  *  1(1=1, 2)  will  be  treated. 

A  more  general  case  can  be  obtained  in  a  similar  manner.  From  (4.15),  let 

the  test  statistic  be 

R  =  Rj  +  R2  +  R3  +  2pxp2[l+  ^  (rj+r2)]  (5.5) 

Rj  =  (m+r2)trB11(B11+S11  )_1  with  =  w2iW21 

-  -1  1 
R2  =  (m+r1)trB22(B22+S22)  with  B22  =  W12W12 


where  m  =  n-r 
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R3  =  m  (m+'"1)(ni+r2)tr(B11+S11)  ls12^ B22+S22^  ls21 


m  (m+r2)trSi]fBii+Sli)  "  m  (m+r1)trS22(B22+S22)  . 

Theorem  5.1.  For  Problem  [II]  :©12  3  0,  e21  =  0,  and  z 12  =  0,  the 
asymptotic  null  distribution  of  R  in  (5.5)  is  given  by 


P(R  1  x)  =  Gf(x)  -  (flSl+f2s2+f3s3)[Gf(x)  -  2Gf+2(x)  +  Gf+4(x)] 


-  ^  P1P2(r1+r2)[2Gf(x)  -  3Gf+2(x)  +  Gf+4(x)]  +  0(m"3/2) 


(5.6) 


where  f  =  f^  +  f2  +  f3>  f^  =  *3lr2*  ^2  =  '32rl’  ^3  =  ^1^2*  S1  =  Pj  +  r2  +  1 

S2  =  p2  +  rl  +  1  and  s3  =  P1  +  p2  +  ** 

The  case  of  Problem  [I]  follows  directly  from  this  theorem. 

Corollary  5.1.  For  Problem  [I]:  e12  =  0  and  e^2  =  0,  the  asymptotic  null 
distribution  of 


T  3  T(r1/P1)  +  +  "2m  rl] 


(5.7) 


with  T(B)  in  (3.22)  is  given  by  (5.6)  with  r2  =  0. 

Outline  of  the  Proof  of  Theorem  5.1.  Under  the  null  hypothesis,  assuming 
z  =  I  without  loss  of  generality,  we  have  the  three  independent  Wishart  variates 
B. ,  -  W  (I,r0),  B00  ~  W  (I,r,)  and  S  ~  W  (I,m).  As  usual  let 

p^  i.  Li.  p2  i  p 


V  =  Vn  Vl2  =  ^  (I  S-I ) 
V21  V22 


(5.8) 


and  expand  R. 's  in  terms  of  V,  B^  and  B22  as 
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"J“fBjl+±,J1+Iq.2  +  Op(m-3''2)(j=l,2) 

Rj  =  trV12V21  -  2PjP2  +  i  q31  +  i  tq32  +  q32  -  PjP2(iyr2)]  +  Op(nf  3/2) 

where  qjl  =  -  trBjjVjj,  qj2  =  trB^V2.  -  trB2j  +  FjtrB^  with  Fj  =  r2  and 

r2  °  rl  <*31  °  "  trVllV12V21  "  trV12V22V21’  q32  =  trVllV12V21  + 

trV12V22V21  +  trVllV12V22V21  and  q32  *  ^rl+r2*trV12V21  *  trBllV12V21  " 

trB22V21V12  +  p2trBll  +  P2trB22"  Then  R  1S  exPressed  as 

R  =  trBn  +  trB22  +  trVJ2V21  +  ~  (pn+p21+q31^ 

(5.9) 

+  i  (q12+q22+q32+q32}  +  2plp2  [1+  S  {rl+r2)]  +  Op(">"3/2)  . 

Then  the  characteristic  function  C(t)  of  R  is  expanded  as 

C(t)  =  E{H(t)  [1  +  J- Ax  A2]}  +0(m"3/2).  (5.10) 

»^m 

where  H(t)  *  exp[it  trBjj  +  it  trB22  +  it  and  ^1  and  ^2  are  ^unct^ons 

of  B^'s  and  V^.'s.  Since  V^.'s  are  not  independent,  in  evaluation  of  the 

expectation  in  (5. 10), we  may  use  the  following  lemma. 

Lemma  5.1.  (1)  The  pdf  of  V  is  expanded  as 

f(V)  =  c  exp(-  j  trV2)(l+  —  [-  |  (p+l)trV  +  |  trV3]}  +  0(m_1). 

i/m 

(2)  The  conditional  pdf  of  V12  given  Vjj  and  V22  is  expanded  as 

f(V12> V11*V22)  =  C  exp("  1  trVl2V21){l+  l~  P2trVll  "  PltrVll 

Zvm 

+  trVnV12V21  +  trV12V22V21]}  +  0(m'1). 


Proof.  (1)  is  well  known.  (2)  is  obtained  by  f(V)/f  (V^)f2(V22) »  where 
the  marginal  pdf's  of  V..'s  are  first  expanded  as  in  (1). 

Lemma  5.2.  The  characteristic  function  of  R  is  evaluated  as 


C(t)  =  (l-2it)*f/2{l  -  ^  (f1s1+f2s2+f3s3) [(l-2it)_1-l]2 

‘  i  plp2(rl+r2)[2  -  3(l-2it)-1  +  (l-2it)'2]}  +  0(nf3/2). 


(5.11 


Proof.  The  proof  is  straightforward  although  it  involves  a  lot  of 
computation.  The  result  is  obtained  by  using  Lemma  5.1  and  the  following  well 

known  results: 

2 

E[exp(1ttrV12V21HH  ^  ,31  ♦  i  (1t)q32  *  q^)» 

*  E  exp(itm  tr  SjJs^S^S^)  +  0(»"'3/2) 

-f  12 

=  (l-2it)  3  [1  -  f3S3{(l-2it)‘1  -  I)2  (5.12) 

+  0(m"3/2)] 

(Fujikoshi  (1970),  Muirhead  (1970)) 


E[exp{it(m+r2)trB11(B11+S11) 


-1 


-1. 


(5.13) 


+  it(m+r1)trB22(B22+S22)"  >] 


(Fujikoshi  (1970),  Murihead  (1970)) 

Inverting  C(t)  in  Lemma  5.2,  we  obtain  the  result  in  Theorem  5.1. 
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6.  TESTS  FOR  INDEPENDENCE  WHEN  Bn  =  0  AND  6^  =  0 


As  has  been  stated  in  Section  2,  for  problem  [III]  H:z12  =  0  given 
9j2  =  0  ancl  02i  =  0  KFK  (1984)  proposed  the  LBI  test,  a  LRT-like  test  and  a 
trace  test,  and  considered  the  asymptotic  null  and  nonnull  distributions  of 
these  tests.  In  this  section,  adopting  what  we  call  the  method  of  small-small 
asymptotics  due  to  Mukerjee  and  Chandra  (1984),  we  compare  the  power  functions 
of  those  tests.  Since  the  model  in  Problem  [III]  is 


(6.1) 

(6.2) 


it  contains  the  MANOVA  model  as  a  special  case  with  Xj  *  X2.  The  comparison  we 
make  here  deals  with  the  comparison  of  the  Pi 11 ai  test  and  the  LRT  for  independ¬ 
ence  in  the  MANOVA  model.  In  fact,  as  will  be  shown,  when  X^  *  X2,  the  LBI 
test  and  the  trace  test  in  Problem  [III]  are  both  reduced  to  the  Pi  1 1 ai  test  of 
independence  in  the  MANOVA  model,  which  is  LBI,  while  the  LRT-like  test  is 
reduced  to  the  LRT  of  independence  in  the  MANOVA  model.  Following  KFK  (1984), 
assume  Pj  >_  p2  without  loss  of  generality  and  let 


Q0  *  I  “  X(X'X)Y  =  lQl'0  with  ZqZq  =  inQ 
Qi  *  I  -  Xi(X*Xi)”1Xi  (1-1,2) 


(6.3) 


where  Zginxng,  n^  =  n-rg  and  rg  =  rank  [X^.X^and  A+  is  Penrose  inverse  of  A. 


Further  let 


x(x'x)V  -  x1(x^x1)"1  -  with  z!z.  =  ir^_r_ 
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where  Z. :nx(rQ-r. ) ,  and  let 


Zi  =  [Zi *Zo] :  nx{n0+r0‘ri)* 


(6.5) 


Mi  (r0‘r1)  =  Z‘Y  =  Zi  Yi 
^iJ(n-r0)  1  1  z'  Y. 


(6.6) 


S  =  G  +  B  and  R  =  S12S22S21Sjj 


(6.7) 


I 

where  S  =  (S..)  with  $..:p.xp  G  *  (G,.)  with  G,.  =  U.U.rp.xp.  and 
IJ  ■  J  ■  J  *  J  1J  1  J  I  J 


B  =  (B..)  for  i,  j=l,2  with 

*  J 


B11  S12 


?21  B22 


HjMj  MjKM^ 

M2K'M1  M2M2, 


,  K  =  Z,Z2 


(6.8) 


Here  we  note  that  when  Xj  *  X2>  Qq  =  Qi  (i=l,2)  so  that  B  *  0  and  S  ~  G 
with  G  '  W(z,ng),  which  is  nothing  but  the  canonical  form  for  the  problem 
of  testing  =  0  ln  the  MANOVA  model.  Now  based  on  the  notation  introduced 
above,  our  problem  is  to  test  Ej2  *=  0  based  on  S  in  (5.7)  and  then  the  LRT- 
like  test  statistic,  the  tract  test  statistic  and  the  LBI  test  statistic 
considered  in  KFK  are  respectively  expressed  as 


Tj  =  -  nQ  log  1 1  -  R| 


(6.9) 


ng  tr  R 


n0  {nin2trR  "  nlp2trSllYlQlQ2C)lYl 
’  n2PltrS22Y2Q2QlQ2Y2>  +  P1P2^2“  n7 


(6.10) 


(6.11) 


where 


ru  =  n  -  r.  (i=0,l,2)  and  t  =  -  2r. 


(6.12) 


As  has  been  shown  in  KFK,  the  power  functions  of  these  invariant  tests 


depend  on  (0,1)  only  through  the  canonical  correlations  pf  >  ...  >  p^ 

1  ~  ”  P2  or 

z  or  the  characteristic  roots  of  2i2s22s21sll *  ^ere  to  cons^er  t^ie  sma^“ 


small  asymptotics  for  the  power  functions,  first  fix  _> 
a),  t  0,  which  are  chosen  to  be  small  later,  and  take 


>  m  >0  with 

-  P2~ 


12-1  (i-1  n  \ 

2  pj  '  nn  “j  (j_1 . p2 


(6.13) 


1  ^ 
where  —  is  supposedly  small.  This  implies  p.  =  (2u./nn)2  is  eventually 

ng  J  J  U 

small -small.  Further  let 


5  =  ay,  +  ...  +  a)  *  tr$  (6.14) 

1  P2 

*  =  di ag-Cm , , . . .  ,wn  1  (6.15) 

1  p2 

f  =  pjP2  and  s  =  Pj^  +  p2  +  1.  (6.16) 


Then  from  the  results  in  KFK,  the  asymptotic  power  functions  of  the  tests 
based  on  T. 's  in  (5.9),  (5.10)  and  (5.11)  under  the  local  alternative  (p  . > 

•  sJ 

in  (5.13)  (as  nQ-»®)  are  given  by 


P(T.  >  x.)  -  Sf(x.6)  +  ^  Zj=Q  aijBf+2j(x1:6)  +0(n"2),  (6.17) 


(i=l,2,3)  where  5k(x:6)  =  1  -  Gk(x:<5),  Gk(x:5)  is  the  distribution  function 
2 

of  X  distribution  with  degrees  of  freedom  k  and  noncentrality  parameter  6, 


38 


a1Q  =-i  fs  -  \  f  (t+trKK* )  -  (trKK*  )<5  -  tr$2 

all  =  7  fs  +  7  f(t+trKK')  "  tfi 

a12  =  (t+trKK  )5  +  2tr$2,  =  ■  tr$2,  =  0 

a20  =  '  ?  fs  "  7  f( t+trKK')  -  (trKK* )S  -  tr$2 

a21  =  7  fs  +  7  f(t+trKK')  "  tfi 

a22  =  "  7  fs  +  (s+t+trKK  )<$  +  2tr$2 

2 

a23  =  -  sS,  a24  =  -  tr$ 

a30  =  "  7  fs  ‘  7  ftrKK'  "  (trKK ' )6  -  tr$2 

a31  =  7  fs  +  7  ftrKK' 

a32  =  ”  fs  +  (s+trKK  )&  +  2tr$2 

2 

a33  =  ~  s6  and  a34  =  "  tr$  . 

Note  that  E^=q  a..  =0  (i=l,2,3).  To  evaluate  the  power  functions  in  (6.17) 
further,  let  g^( x : 6 )  be  the  pdf  of  G^(x:6)  and  let  Gf(x)  *  Gf(x:0)  and 
gf(x)  *  gf(x:0).  Then  it  is  easy  to  see  that 

Gf+2(x:5)  =  2qf+2(x:5)  +  Gf(x:<5)  (6.18) 

gf+2(x)  =  xgf(x)/f  (6.19) 

5f(x+  jjj:6)  =  Gf(x:6)  -  £  gf(x:6)  +  o(nf2).  (6.20) 

Using  (6.18),  the  power  functions  are  evaluated  as 


where 


P(Ti  >  x.)  =  £f (x. :fi)  +  ^  Ej=1  bij9f+2j(xi :6)  +  °(n02)  (6*21) 

bll  =  \  fs  +  fU+trKK*)  +  2 ( trKK ' )s  +  2tr$2 
b12  =  2(t+trKK* )6  +  2tr$2 
-  -  2tr$  t  b^4  =  ® 

b2i  *  \  f S  +  f(t+trKK* )  +  2(trKK')5  +  2tr$2 
b22  =  -  |  fs  +  2(t+trKK* )<S  +  2tr*2 
b23  =  -  2s6  -  2tr$2,  b24  =  -  2tr$2 

b31  *  \  fs  +  f  trK1C'  +  2(trKK')«  +  2tr$2 

b32  =  "  1  fS  +  2(trKK')6  +  2tr®2 
b33  =  -  2s6  -  2tr$2 ,  b34  =  -  2tr$2. 

Using  (6.21)  and  (6.19),  under  the  null  hypothesis  6=0 

a  =  P(T.  >  x.|H0)  =  Gf(x.)  -  ±  bi0gf(x.)  +  0(n'2)  (6.22) 

(i=l,2,3)  where 

b10  *~*{l  s+t+trKK')  =-  x1C1 

sx 

b20  =  "  x2(i  s+t+trKK  •  ?(y+l ) }  -  ■  X2C2  (6.23) 

~  1  '  sx3 

b30  =  "  X3^I  s+trKK  "  2( f+2) ^  -  “  X3C3 


From  (6.20)(or  Hill  and  Davis  formula)  we  obtain  with  u  in  G^(u)  =  o 


xi  =  u  “  nQ  bi0  +  °'n02^ 


(6.24) 


where  =  -  uC..  with  in  (6.23).  Hence  using  (6.20)  and  (6.24),  we 
obtain  the  following  theorem. 

Theorem  6.1.  Let  ir(<j>..  ,$)  be  the  power  function  of  the  test  <|>.  of  size 
a  with  critical  region  >  x.(i=l,2,3).  Then  for  u  satisfying  G^(u)  =  a 
and  for  6  small,  it  is  evaluated  as 

ir(*i.*)  -  Gf(u:6)  +  ^  H.(s)gf(u)  +  0(n‘2)(i=l,2,3)  (6.25) 

where 


H^S)  =c.j6  +  ci262  +  c13tr$2  +  0(nQ2)(i=l,2,3) 


(6.26) 


cn  ■  ca  '  c3i  -  T  trKK'  -  nky 


c12  =  '  T  trKK'  +  (s+2trKK’>  - 


C22  =  C32  =  °12  +  f(f+2)2(f+4)(f+6) 


.  _  2u  .  2uc 

c13  "  T  +  7T7+2T 


C23  =  C33  =  C13 


Proof.  Using  (6.20)  and  (6.24),  the  power  functions  in  (6.21)  are 
directly  shown  to  be  equal  to  those  in  (6.25)  after  some  algebra,  where 
gf+3j(u:6)  -  e'5  r^=0(6k/kI )gf+2j+2k(u)  and  gf+2(u)  =  ugf(u)/f  were  used. 
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Pp 

In  the  expression  (6.25),  6  has  been  chosen  to  be  small.  That 

is,  the  power  functions  were  first  asymptotically  expanded  in  the  orders  of 

-k  -1 

Oq  (k=0,l,2...),  and  then  for  6  small  the  terms  of  order  n^  were  asymptot¬ 
ic  0 

ically  expanded  in  the  orders  of  6  (k=0,l,...)  since  the  terms  of  order  ng 

are  common  to  all  the  tests,  i.e.,  S^(u:6).  From  this  expression,  the  local 

behaviors  of  the  power  functions  are  compared  as  follows. 

Theorem  6.2.  For  the  tests  <^ ,  it  follows  that  with  6  =  tr<t 

(1)  lim  lim  nn[7rU . ,4>)  -  *(♦.,«)]/«  =  0  (j=2,3) 

6+0  n0*»  u  J  1 

(2)  lim  lim  nn [ir (<#» . ,$)  -  Tr(^1  =  a(4»)(j=2,3) 

6+0  nQ+*>  u  J  1 

and 

(3)  lim  lim  nQ [tt(4>3 =  0 
6+0  n^-*® 

where  a($)  =  [2u^g^.(u)/f(f+2)^(f+4)(f+6)]Y($)  with 

y(*)  =  (p,+p?+l)  -  (p,p?+2)  lim  [tr$2/(tr$)2]  (6.27) 

1  c  1  ^  6+0 

Proof.  Immediate  from  Theorem  6.1. 

Now  (1)  implies  that  in  terms  of  power  all  the  three  tests  are  asymptotically 
equivalent  up  to  o(ng*)  and  0(6  }.  The  asymptotic  difference  between  the  LBI 

test  <|>3  (or  the  trace  test  $2)  and  the  LRT-like  test  <j>j  appears  in  the  term 

-12  2  2 
of  o(nn  )  and  0(6  )  as  is  shown  in  (2).  Setting  x  =  lim[tr$  /(tr$)  ],  from 

u  6+0 

(2),  if  y($)  >  0,  or  equivalently  (p1+p2+l)/(p^P2+2)  >  x,  the  LBI  test  is 

-1  2 

asymptotically  better  up  to  o(ng  )  and  o(<5  )  than  the  LRT-like  test,  while  if 

2  2 

y(*)  <  0,  the  LRT-like  test  is  asymptotically  better.  Since  tr$  =  Zu^.  and 
tr$  =  Zoo.  =  6,  the  inequality 

Zo)2  <  [Zco.]2  <  p2[Eu?] 


follows  from  w.  >  0  and  Schwartz's  inequality.  This  implies 

-i-  <  t  =  lim  tr$2/[tr$]2  <  1  (6.28) 

p2  "  5-0 

The  equality  in  the  first  of  (6.28)  holds  if  and  only  if  =  to 2  -  . . .  =  o> p  , 

while  the  equality  in  the  second  inequality  holds  if  and  only  if  u*  =  ...  * 

oj  =0  since  Ui  >  w  o  >  • . .  >  u  „  >0.  On  the  other  hand, 

P2  1  —  2  —  —  p2  — 

~  <  Tq  =  (P1+P2+1)/(P1P2+1)  <  1 

since  Pj  1  P2  51  0-  Hence  both  the  cases  1  t  >_  tg  and  tq  >  t  >_  l/p2  can 

occur.  The  above  observation  will  show  that  the  closer  t  is  to  1,  the  more 


concentrated  m. 1 s  are  around  ^  =  ... 
the  more  spread  w. 's  are. 


to  ,  while  the  closer  t  is  to  l/p?, 
P2 


Next,  we  consider  the  small-small  asymptotic  comparison  between  the  LRT 


and  the  Pi 11 ai  test  of  Independence  in  the  MANOVA  model.  As  has  been  pointed 
out,  the  problem  of  testing  independence  in  the  MANOVA  model  is  included  as 
a  special  case  of  our  problem  with  =  Xg.  In  case  X^  =  X2,  =  0  in  (6.4) 

so  that  S  =  G  in  (6.7),  K  =  0  in  (6.8),  nj  =  n2  =  no  in  (6*12),  and  T2  =  Tg 
in  (6.9)  and  (6.11).  However  this  does  not  cause  any  changes  in  the  results  of 
Theorem  6.1  and  6.2  except  for  the  slight  changes  of  the  coefficients  c..'s  in 

■  J 

Theorem  6.1.  That  is,  by  setting  K  =  0  in  Cjj  and  c12»  results  in  Theorem 

6.1  holds  as  it  is,  while  Theorem  6.1  is  effective  whether  or  not  Xj  =  X^. 

Corollary  6.1.  For  testing  independence  in  the  MANOVA  model  with  X1  -  X2, 

all  the  results  in  Theorem  6.1  hold.  If  Tq  >  t(>^  I/P2) *  Pillai's  test,  which 

-1  2 

is  LBI  for  fixed  ng,  is  asymptotically  better  up  to  0(n“  )  and  0(5  )  than  the 
LRT  while  if  1  >  t  >  ig,  the  LRT  is  asymptotically  better. 
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7.  REMARKS 


In  this  section,  we  first  consider  Problem  [IV]:e12  =  0  and  e12  ■  0. 

For  this  problem,  it  is  easy  to  see  that  a  canonical  form  of  the  model  is 

also  given  by  the  model  in  (4.2)  where  the  joint  hypothesis  e12  =  0,  e21  =  0 

and  e,2  =  0  in  Problem  [II]  was  tested,  and  the  hypothesis  here  becomes  n12  ~  0 

rlpl  r2p2 

and  o21  =  0.  Further,  the  same  group  <S  =  (^(p^xG^p^xR  xR  acting  on 
(W,S)  as  in  (4.5)  leaves  this  problem.  Hence  the  class  of  invariant  tests  in 
Problem  [II]  is  exactly  the  same  as  the  class  of  invariant  test  in  Problem  [IV] 
we  are  presently  considering.  This  implies  that  the  power  function  of  an 
invariant  test  in  Problem  [IV]  is  also  a  function  of  n,  £j2  and  ?21  in 
However,  under  the  null  hypothesis  C^2  =  0  and  =  0,  it  still  depends  on 

n  because  n  may  not  be  zero  in  Problem  [IV].  Therefore  in  general  an 
invariant  test  in  the  present  problem  is  not  similar.  In  fact,  this  follows 
from  the  fact  that  the  group  does  not  act  transitively  on  the  parameter  space 

I 

of  the  null  hypothesis.  This  is  true  even  in  the  case  X^X2  =  0  where  the 
group  is  enlarged  to  3  *  ^(r^x 0(r2)xG  as  in  (  ).  For  example,  suppose  we 

construct  such  statistics  as 


^2  =  l'*22^'*22  +  ^12^12^  and  ^"3  =  ^11^^11  +  ^21^21^  (^*1) 

for  n12  =  0  and  n21  =  0  respectively  analogously  to  (4.8),  or 
L2  =  trW12S22W12  and  L3  =  trW21SllW21' 


But  here  L2  and  (or  L2  and  L^)  are  correlated  under  the  mull  hypothesis  so 

I  t 

that  any  test  combining  L2  and  (or  L2  and  L^)  is  not  similar  unless  one  of 
the  two  statistics  is  completely  ignored.  Because  of  the  non-similarity 


feature  of  the  problem  we  leave  it  here.  One  might  use  a  non-similar  test 
by  combining  L 2  and  L3  in  (7.1)  in  such  a  way  as  L2L3  or  might  test  the  two 
hypothesis  separately.  It  is  noted  that  an  explicit  form  of  the  LRT  for  the 
present  problem  is  difficult  to  derive. 
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